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INTRODUCTION 



, Rotating relativistic stars are among the most promis- 
5_j ' ing sources of gravitational waves (see Sathyaprakash and 
Schutz 2009 for a review), especially when they are newly 
born, be it as proto-neutron stars resulting from a gravita- 
tional core collapse or as compact remnants of neutron stars 
binary mergers. Yet, their study is a complicated task, in- 
volving the solution of Einstein's equations in a dynamical 
regime while, at the same time, many microphysical phe- 
nomena need to be taken into account. Therefore, one of 
the first compulsory steps in their modelling is the analysis 
of axisymmetric and stationary configurations of rotating 
self-gravitating perfect fluids. However, even when a con- 
stant rotation profile is assumed, the solution space of this 
academic problem is known to possess a complex structure 
(see Ansorg et al. 2004; Meinel et al. 2008), which implies 
that robust and precise numerical codes are needed before 
one can proceed with the astrophysical investigation. 

A well-studied example is the solution space of homo- 
geneous, uniformly rotating bodies to which, in Newtonian 
gravity, the famous Maclaurin spheroids belong. For given 
mass density, this family depends, but for a scaling factor, 
on only one parameter, for which reason we may refer to it 
as the Maclaurin sequence. However, the full solution space 
is much richer as it contains infinitely many Newtonian se- 
quences that branch off from the Maclaurin sequence (see 



ABSTRACT 

A highly accurate, multi-domain spectral code is used in order to construct se- 
quences of general relativistic, differentially rotating neutron stars in axisymmetry and 
stationarity. For bodies with a spheroidal topology and a homogeneous or an A'^ = 1 
polytropic equation of state, we investigate the solution space corresponding to broad 
ranges of degree of differential rotation and stellar densities. In particular, starting 
from static and spherical configurations, we analyse the changes of the corresponding 
surface shapes as the rate of rotation is increased. For a sufficiently weak degree of 
differential rotation, the sequences terminate at a mass-shedding limit, while for mod- 
erate and strong rates of differential rotation, they exhibit a continuous parametric 
transition to a regime of toroidal fluid bodies. In this article, we concentrate on the 
appearance of this transition, analyse in detail its occurrence and show its relevance 
for the calculation of astrophysical sequences. Moreover, we find that the solution 
space contains various types of spheroidal configurations, which were not considered 
in previous work, mainly due to numerical limitations. 

Ansorg et al. ^OoM))- 

Taking these configurations into the 
realm of general relativity, one recognizes the entire com- 
plexity of the solution space which turns out to consist of 
infinit ely many disjoint solution classes (see Ansorg et al. 
(j2004l )). Within these classes (which are two-dimensional 



subsets of the solution space) , specific sections of the above 
Newtonian sequences form particular limiting curves. 

One consequence of this class structure is that any 
continuous sequence of relativistic solutions starting at a 
static and spherical body and running through the asso- 
ciated spheroidal class, is not connected to the regime of 
configurations with toroidal shape. Instead, such a sequence 
eventually leads either to mass shedding limi10 or to a limit 
of infinite central pressure. 

In a slightly more realistic model of a star, which one 
obtains by taking polytropic matter, a picture similar to the 
one de picted above for homogeneous bodies emerges (Meinel 
et al. (|2008l )). As one considers polytropes with increasing 
polytropic index A'' (with A'' = describing the homogeneous 
bodies), one finds that the different classes drift apart from 



^ A mass-shedding limit is given if a fiuid particle at the surfaee of 
the body moves with the same angular veloeity as a test particle 
at that spatial point. The corresponding geometrical shape of the 
fiuid body possesses a cusp there. 
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one another, implying an even clearer separation between 
them. 

However, apart from the equation of state (EOS), actual 
stars differ from the model of simple homogeneous rigidly ro- 
tating fluids through their presumable differential rotation, 
which is expected to occur in proto-neutron stars born from 
a gravitational collapse or from the merger of binary neu- 
tron stars. Hence, a natural question that arises is that of the 
influence of the rotation profile on the main characteristics 
of a stationary, relativistic star, such as its maximal mass 
and angular momentum, but also its actual geometrical sur- 
face shape. The evolution of the latter one possibly leads to 
various types of instabilities and corresponding signals (e.g. 
electromagnetic or gravitational waves, neutrinos). Due to 
its importance for estimating the delay between a neutron 
stars binary merger and the collapse to a black hole, this 
issue was already the topic of various studie s, be i t for cold 
or warm equations of state (see Stergioulas (|2003l ) for a re- 
view on rotating stars in relativity). However, in the case 
of differentially rotating stars, most of the work done up to 
now seemed to suffer from numerical limitations when the 
rotation profile started to be too sharp. Moreover, none of 
those studies really dealt with the changes undergone by the 
solution space due to the rotation law. 

In this article, we focus on the influence of differen- 
tial rotation on the parameter space of rotating relativis- 
tic star^. In particular, we only consider star-like bodies, 
i.e. configurations with spheroidal topology (simply con- 
nected ones), and use two basic EOSs (homogeneous and 
polytropes with = 1) as well as one specific type of ro- 
tation law. The latter one is the most classical in general 
relativity, introduced by Komatsu et al. (|l989l ). in which 
a parameter A appears that measures the degree by which 
the rotation is differential. This allows us to go from rigid 
rotation to configurations with a strong gradient in the ro- 
tation profile. As will be described below, in this situation 
some critical values of A arise, at which a merging of typical 
classes is exhibited. This result shows that the introduction 
of several disjoint classes, initially done in order to charac- 
terize the solution space corresponding to uniform rotation, 
only makes sense for sufficiently weak differential rotatiorjfl. 

The paper is organized as follows. First, we clarify the 
distinction between the so-called spheroidal class and the 
toroidal class of rotating stars, a notation that should not 
be confused with the star's topology, see discussion below. 
Then, in Section 3, we introduce the Lewis-Papapetrou line 
element that describes stationary and axisymmetric space- 
times. At the end of that Section we discuss the Komatsu law 
of differential rotation, to be used in this paper, as well as 
the resulting general relativistic field equations. In Section 4 
we briefly describe particular features of our pseudo-spectral 
code, used to obtain highly accurate numerical solutions and 

^ We postpone the discussion of astrophysically relevant quan- 
tities (such as the mass, the angular momentum) to some other 
article. 

^ Hereby, these critical values of A depend sensitively on the EOS 
being chosen. Below we demonstrate that for N = 1 polytropes 
the class structure can still be introduced for sufficiently weak dif- 
ferential rotation. In contrast, for homogeneous matter the class 
structure ceases to exist in a strict sense for any rate of differential 
rotation. 



present numerical results. We introduce four speciflc types 
of sequences of differentially rotating relativistic stars, which 
describe the solution space in question independently of the 
EOS being chosen. Finally Section 5 provides a summary 
and a discussion of our results. 



2 SPHEROIDAL AND TOROIDAL CLASSES 

In the following we shall concentrate on two specific classes, 
which we call spheroidal and toroidal. For rigidly rotating 
stars, the first one has already been mentioned in the Intro- 
duction and contains continuous parametric sequences that 
start at a static and spherical body. The shapes of these 
bodies possess spheroidal topology throughout the entire 
class. Any such sequence eventually leads either to a mass- 
shedding limit or to a limit of infinite central pressure. For 
this reason it makes sense to identify extremal values for cer- 
tain physical quantities, such as the maximal baryon mass, 
within the spheroidal class. 

For the introduction of the toroidal class we consider 
the so-called Dyson ring sequence, which is the Newtonian 
sequence of rigidly rotating homogeneous bodies branching 
off from the Maclaurin sequence at the branch point, at 
which the Maclaurin spheroids become secularly unstable 
with respe ct to the first ax isymm etric perturbation (Chan- 
drasekhar (|l967l ). Bardeen (|l97ll )). One branch of the Dyson 
ring sequence terminates at a mass-shedding limit, while an- 
other one leads to the regime of figures with toroidal topol- 
ogy and finally to arbitrarily thin ringfl. Taking the Dyson 
ring solutions into the relativistic regime and choosing again 
homogeneous matter, we find an associated toroidal class 
of solutions which is disjoint with respect to the spheroidal 
class. In contrast to the spheroidal class, it possesses contin- 
uous sequences of relativistic solutions leading to infinitely 
thin tori. Moreover, the solutions of the toroidal class do not 
possess maximal possible values of the baryon and gravita- 
tional masses. Instead, those quantities increase to arbitrar- 
ily high values as the relativistic tori become thin ner an d 
thinner (see Ansorg et al. (|2003d ) and Fischer et al. (|2005h ). 

As explained in the Introduction, the solution space of 
rigidly rotating A'^ = 1 polytropes is quite similar, except 
that even in the Newtonian limit there is already a clear 
gap between the spheroidal and toroidal sequences. This gap 
grows further as polytropes with increasing index N are con- 
sidered. 

As an important remark, note that in the above def- 
inition of spheroidal and toroidal classes, the existence of 
continuous transitions, either to spheres or to infinitely thin 
tori, is considered. This does not determine a priori the ge- 
ometrical topology of the configurations which form the se- 
quence. This is to say that a body with spheroidal topology 
(i.e. which is simply connected and has no hole) may in fact 
belong to the toroidal class, provided that a continuous pa- 
rameter transition to the regime of infinitely thin tori can 
be identified. 

In this article we focus, for the purpose of simplicity, on 
stars with spheroidal topology and find that the situation 

^ We call a ring thin if the ratio of inner to outer circumferential 
equatorial radius is close to 1. 
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becomes more complicated for differentially rotating bod- 
ies. In particular, a merging of the spheroidal and toroidal 
classes is exhibited, meaning that from some critical rate 
of differential rotation on the classification does not make 
sens43- However, as we shall see below, a discussion of the 
corresponding solution space is possible in terms of four dif- 
ferent types (A,B,C and D) of one-dimensional parameter 
sequences. 



3 DIFFERENTIALLY ROTATING 

SPHEROIDAL STARS IN GENERAL 
RELATIVITY 

3.1 Line element 

Using Lewis-Papapetrou coordinates {Q,z,ip,t), the line el- 
ement for a stationary, axially symmetric spacetime corre- 
sponding to a rotating perfect fluid conflguration can be 
written in the form (see e.g. Stephani et al. ()2003l )^ 



As' 



e''''(dfj" 4- dz") + W^^e-^"(d<y3 - lj Atf - e^" di", (1) 



where the metric functions /x, W , v and u) are functions of 
Q and z alone. Note that 



VK = on the rotation axis q = 



(2) 



The coordinates (g, t) are uniquely defined by the re- 
quirement that the metric coefficients and their first deriva- 
tives be continuous everywhere, in particular at the surface 
of the rotating body. 



3.2 Equation of state and Rotation law 

For a barotropic perfect fiuid, the energy-momentum tensor 
reads 



(3) 



where is the fiuid's 4-velocity. The total mass-energy den- 
sity e as a function of the pressure p is given via an equation 
of state. 

In this paper we consider neutron stars described by 
polytropic or homogeneous EOSs. The polytropic equation 
of state is given by 



e{p) =p + \JpIK 



(4) 



where K is the polytropic constant. The polytropic index TV 
is chosen to be unity as this value fits well most of th e real is- 
tic EOSs for neutron stars (see Lattimer & Prakash (|2006l )'). 
Since this is no longer true for strange stars (e.g. Limousin 
et al. (|2005l )). whose properties can be much better repro- 
duc ed us ing a homogeneous mo del (se e Gondek-Rosihska et 
al. (|2003l '). Amsterdamski et al. (|2002l )'). we also use such an 
EOS with 



e{p) — eo — constant. 



(5) 



^ Note that the numerical code needs to be specifically adapted 
to the treatment of configurations with rotation rates close to a 
critical value, in order to avoid non-convergence or random nu- 
merical jumps between different types of solutions. 



The stationary and circular motion of the fluid particles 
implies a 4-velocity of the form: 

u' ^e-'^ie +iiv') (6) 

where 5' and 77' are the Killing vectors with respect to 
stationarity and axisymmetry respectively. The coefficient 
= u'/ti*' denotes the angular velocity of the fluid and 
depends on q and z. The integrability conditions of the Ein- 
stein equations, i.e. the Euler equations 



r'^fc = 



(7) 



restrict the possible form that the angular velocity can take 
They imply 

dp 



e+p 



dV + u u,j,AQ. 



(8) 



from which it follows that the product u is a function 
depending on SI only. 



and that the specific enthalpy 



h = 



e + p 



(where es denotes the baryon mass density) with 

dh dp 
h e + p 

is related to the metric quantities as follows: 

h = ho exp ^o-V - J F{Q,)dflj . 



(9) 
(10) 

(11) 
(12) 



Here, the constant Sic denotes the central angular velocity, 
while Vo, i.e. the value of the function V at the fluid's north 
pole {g ~ 0, z = Vp), is a. relativistic parameter charac- 
terizing the rotating fluid body. It is related to the relative 
redshift zq of photons emitted from the north pole and re- 
ceived at infinity via 

-Vo 1 

zo =e ° - 1. 

Note that from u^Ui = — 1 we have 
1 
2 

and 

/2 



(13) 



w'{n-uj) 



(14) 



Hence, by virtue of ((Qjl we obtain Q, in terms of the metric 
coefficients as the solution of the implicit equation 



41. _ M^2(J^ _^)2- 



(15) 



It follows from Q and jlS} that (i) F{Q.c) = and (u) 
SI = S7c on the rotation axis. 

The enthalpy ho at vanishing pressure takes, for ordi- 
nary matter (such as the ones considered here), the value 

^ An exceptio n is so called strange quark matter, see e.g. Gour- 
goulhon et al. lll999l') . 
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The function F can be chosen freely, being only subject 
to the condition -F(Slc) ~ 0. In this paper we use the simplest 
form which goes beyond uniform rotation, i.e. we take a 
linear ansatz 



F{n) = (fie - n) , 



(16) 



with the constant A describing the rate of differential rota- 
tion. This particular rotation law, originally introduced by 
Komatsu et al. (|l989l ). has been considered by vari ous au - 
thors (including Cook et al. (| 19921 ). Bonaz zola e t al. (|l993h . 
Go ussard et al. (| 19981 ). B aumga rte et al. (|2000l ). L yford et 
al. (|2003l ). Morrison et al. (|2004l) . and Villain et al. |200i)). 
Since the typical scale of length of the problem is the radius 
of the star and since the rotation profile becomes uniform in 
the limit A oo, we sh all not use A directly, but following 
Baumgarte et al. (|200Cll ) , we parameterize the sequences and 
the degree of differential rotation using the parameter 



A = A- 



Tc/A, 



(17) 



where re is the star's equatorial coordinate radiufl With- 
out ent ering more into the detail of the law (|16|l [see Ster- 
gioulas (|2003l ) for further information], we note that it ver- 
ifies Rayleigh's criterion for st ability against axisymmetric 
perturbations (Komatsu et al. (| 19891 )) and that, in the New- 
tonian limit, it reduces to 



A^p 



(18) 



which can be shown to be in reasonable agreement with 
configurations of newly born neutron stars resul ting from 
core collapse (see Section 4 of Villain et al. (1200-d)) or from 
dynamical merger simulations (see Lyford et al. (|2003l )). 



3.3 Einstein equations and the free boundary 
value problem 

The set of equations, following from Einstein's field equa- 
tions Rik — \Rgik ~ STrTik for the metric in the form ([l]), 
with ((31) and Q, c an be written as follows (see e.g. Bardeen 
& Wagoner (|l97ll )): 

V ■ (BViy) - ifj2B^e-^"'(Vcj)'' = 



I — v' 



(19a) 




I 



0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



Figure 1. Example of a differentially rotating polytropic body 
(with polytropic index A'^ = 1) close to the transition of spheroidal 
to toroidal topology. The parameters of this configurations are 
given hy A = 0.9, log /imax = 0.6 and fp/re = 0.005. In the 
corresponding diagram in fig. [3] this body belongs to the se- 
quence which terminates in a mass-shedding limit (/3 = 0) at 
Tp/re 0.25. The plot shows isocontours of log /i in a meridional 
cross section. The strongly pronounced cusp at the equatorial 
edge indicates that this configuration rotates almost at the mass- 
shedding limit. 



with [cf. (0)] 



B 



W/g and v 



(19e) 



Here the operator V has the same meaning as in a Euclidean 
three-space in which g, z and ip are cylindrical coordinates. 

The equations (|lll - 116|) allow us to express p, e and O. 
and hence all source terms in (|19|l with respect to the metric 
potentials and the three constants fic,^ and Vb- Thus we 
obtain a complete set of elliptic equations for the four metric 
coefficients \i, B,v and uj. 

Note that the fiuid's surface shape is unknown and 
must be determined in the context of the numerical solution 
scheme. Its definition is given through vanishing pressure 
(i.e. through h = ho), which leads by virtue of (|12p to the 
surface condition: 



V = Vo- J F{Q.)dn 

The asymptotic conditions at infinity: 

fj, — u = uj — 0, B = 1, 



(20) 



(21) 



complete the free boundary value problem to be solved. For a 
given EOS, the corresponding solutions depend on the three 
parameters flc,A and Vb- 



V ■ (g^B'^e-^'^S/ij) = -IGTvgB^e^'''^'' (e + p) 



V ■ (gVB) = IGngBe^^p 



ti,eB + l^,a- -^.e + ^^ {^^ ~ B ^VB) 



l-?;2' 

(19b) 

(19c) 



ig2B^e-^''(Va;)2 -47re2^(e-fp), (19d) 



^ Notice that we shall not use the A variable, but mention it so 
as to make a comparison with articles from other authors easier. 



4 NUMERICAL RESULTS 

4.1 Numerical scheme and parameter 
prescriptions 

We use a pseudo-spectral double-domain method to com- 
pute highly accurate numerical solutions corresponding to 
rotating relativistic stars. This techniqu e has been described 
in detail in Ansorg et al. (|2002l : l2003al ). see also Meinel et 
al. (|2008l ). Here we develop the scheme further to include 
the calculation of differentially rotating stars. In particular, 
strongly pinched objects can be computed (for an example 
see fig. [TJ . 

As a special feature, the code permits the prescription 
of three independent parameters which may be chosen for 
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convenience. For example, it is possible to prescribe the pa- 
rameter triple (/imax, r-p/ro, ^) where ftmax denotes the maxi- 
mal enthalpy of the star. Note that this maximal value need 
not be assumed in the coordinate origin q — z — but 
can, for stars sufficiently close to the toroidal regime, be lo- 
cated at {Qh,z = 0) with Qh > 0, see e.g. fig. [T] The dotted 
line on each diagram of figures |3] and [5] separates config- 
urations with the maximal enthalpy located in the stellar 
center (on the right side) from those having ftmax outside 
the coordinate origin. This curve is defined by the vanish- 
ing of the enthalpy's second derivative with respect to q 
at the origin, i.e. through (d^h/dg^)c = 0. We determine 
such a configuration through the prescription of the triple^ 

(/imax, r-p/r-o, (d^/l/de^jc) or (/imax, A, (d^fe/df(^)c). 

Another possible parameter prescription is the triple 
(ftmax, fp/j^e, /?) wh cre the shedding parameter (3 is defined 
as in Ansorg et al. (|2003al '). 



d(p2 



(22) 
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A = 0.4 










A = 0.7 




typo B V 




type C 


A = 0.8 










A = 1.0 












;|| type D i 






/type A 




0.25 




0.5 


0.75 ] 



r-p/re 

Figure 2. Type A-, B-, C- and D-sequences of differentially ro- 
tating polytropic fluid bodies (polytropic index N = I) with pre- 
scribed maximal enthalpy log/imax = 0.2. Shown is the depen- 
dency between the shed parameter /3 and the ratio rp / Tc of polar 
to equatorial coordinate radius for several values of the rotation 
law parameter A. The bold curve corresponds to the scparatrix 
sequence with A = A^rit = 0.75904. 



with z = Zb{Q) describing the star's surface. This parameter 
is chosen such that /3 = in the mass-shedding limit and 
13=1 for (Maclaurin) spheroids. If a sequence exhibits the 
transition to the toroidal regime (i.e. rp/ro — > 0), /3 tends to 
infinity as this transition is encountered. For this reason we 
prefer to consider the rescaled shedding parameter 



(23) 



which tends to in the mass-shedding limit, to i for the 
Maclaurin spheroids as well as in the spherical non-rotating 
limit and to 1 at the transition point to the toroidal regime. 

For the exact calculation of the critical rotation law pa- 
rameter ylcrit, marking the onset of a continuous parametric 
transition of non-rotating spherical stars to toroidal config- 
urations, the parameter triple (/imax, '"p/''o, /?) proves to be 
particularly useful. At fixed value /imax, the parameter A 
can be interpreted as a function of the two variables r^/rc 
and /3. As will be discussed below, this function possesses a 
saddle point with the value A = Acrit- Hence it can be found 
by solving the equations 



dA 



dA 

df3 



(24) 



for the corresponding parameters (rp/rc)crit and /3crit. 
The associated A assumed for the parameter triple 
(ftmax, (rp/re)crit,/3crit) IS the Critical value Acrit in ques- 
tion. In order to determine Aa-it accurately, we compute the 
function A = A{rp/ra, $) in a vicinity of the critical point 
and perform a numerical differentiation with respect to the 
two variables. We then solve the system (|24p by means of a 
Newton- Raphson scheme. Note that again a pseudo-spectral 
treatment ensures a high accuracy of the value Aait obtained 
in this manner. 



4.2 Polytropic bodies 

In order to illustrate the evolution of the solution space as 
the rate A of differential rotation is increased, we consider 
at first fairly weak relativistic polytropic solutions (with 
polytropic index = f) with a fixed maximal enthalpy, 

log /imax = 0.2. 

For rigid rotation, i.e. A = 0, the sequence starting 
from static and spherical bodies with Tp/re = 1 runs through 
the spheroidal class and terminates in a mass-shedding limit 
characterized by /3 = 0, see fig. (2] We assign type A to 
sequences of this kind and note that the spheroidal class is 
made up of such sequences. 

Associated with the sequence of type A with A = 0, 
there is a sequence that contains only configurations with 
toroidal shape and belongs to the toroidal cl ass. I t runs 
within class II in figure 3.29 of Meinel et al. (|2008l 'l from 
the mass-shed curve to the extreme Kerr limit. Since in fig. 
[5] only sequences with spheroidal topology (i.e. rp > 0) are 
displayed, this sequence does not appear here. 

However, as we move to A = 0.4 we now find this 
sequence in the left part of the figure. It is clearly sepa- 
rated from the corresponding sequence of type A which runs 
through the spheroidal class (in the right part of the figure). 
This means that, for A = 0.4, the sequence belonging to 
the toroidal class stretches into the regime of bodies with 
spheroidal topology, i.e. the toroidal class contains config- 
urations with spheroidal topology. We assign type B to se- 
quences of this kind, that appear for A greater than some 
threshold value to be studied in a forthcoming article. In 
the corresponding situation for A = 0.7, the two sequences 
of types A and B still coexist, but they have moved towards 
one another. The merging of the two sequences appears at 



* Note that for triples of prescribed parameters, which do not 
contain the rotation law parameter A, the corresponding value of 
this quantity is determined in the numerical calculation. Likewise 
the parameters Vb and Qc are part of the computational result. 
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the critical value A^it ~ 0.75904. The further increase of 
A leads again to two separate sequences of types C and D. 
The type C sequences show a continuous transition from 
the regime of spherical (with rp/r^ = 1) to that of toroidal 
bodies (at r^/rc = 0). The type D sequences possess mass- 
shedding limits (with /3 = 0) at both ends. As these two 
ends tend towards one another as A is increased, this type 
of sequence exists only for a small range of values of the pa- 
rameter A. This range has the lower bound ^crit and another 
threshold value as an upper bound (to be studied elsewhere) . 
At this upper bound the sequence of type D degenerates to 
a single point describing a particular mass-shedding fluid 
body. By contrast, the type C sequences exist for all values 
A> A 

crit . 

The qualitative behaviour for log /imax = 0.2 is typical 
for the entire solution space. In fig. |3] we show sequences 
for A e {0, 0.1, 0.2, . . . , 1.5} and log Vax £ {0, 0.2, 0.4, 0.6}, 
log/imax = meaning the Newtonian limit. The lowermost 
curves on the right part of the diagrams always correspond 
to rigid rotation A = 0, while the uppermost ones represent 
the configurations with the highest rotation parameter con- 
sidered, A — 1.5. Also shown (in bold) are the separation 
sequences corresponding to the values ^crit. From fig. [2] it 
becomes apparent that at the separation point the function 
A = A{rp/re, (3) has a maximum with respect to the vari- 
able rp/vc and a minimum with respect to the variable /3, 
i.e. a saddle point at which (|24p holds. Moreover the curves 
described by {d^h/dQ^)c = are plot (dotted lines), which 
indicate that, for configurations to the left of this curve, the 
maximal enthalpy is assumed outside the coordinate origin. 

Note that for log /imax = 0.6 two branches of the separa- 
trix sequence (for which Aait = 0.77804) run into the regime 
of toroidal figures, as opposed to only one such branch for 
log /imax G {0,0.2,0.4}. As a consequence, the neighbour- 
ing sequences with shedding ends (i.e. in the diagram for 
log /imax = 0.6, type D curves with A — 0.8 and A = 0.9) 
possess a branch running into the toroidal regime. However, 
the type D sequence for ^4 = 1 is still solely formed of figures 
with spheroidal topology. 

The evolution of the parameter Aait as a function of 
log /imax is displayed in fig. |4] From this figure it can be 
seen that for a sufficiently large rate of differential rota- 
tion, A > A^f^ — 0.71085, the spheroidal and toroidal 
classes merge and hence a continuous parametric transition 
of spherical to toroidal differentially rotating stars is exhib- 
ited. The merging of the classes occurs for A = j4j^j^^^ at 
log /ii^ax = 0.3845. 

4.3 Homogeneous bodies 

For homogeneous and rigidly rotating fluid bodies in Gen- 
eral Relativity, the Newtonian limiting sequences have been 
found to be the separatrix sequences that separate the 
general-relativistic solution classes from one another, see 
Ansorg et al. (j2004l l. figs 4 and 5 therein. Specifically, the 
Maclaurin spheroids with Vp/rc £ [0.17126, 1] form one 
branch of the separatrix dividing the spheroidal from the 
toroidal class. Other branches of this particular separa- 
trix are formed by the Maclaurin spheroids with r-p/rc G 
[0.11160,0.17126] and the (non-relativistic) Dys on rin g se- 
quence (marked by Af in fig. 4 of Ansorg et al. ()2004l ll. 
In the diagram for log /imax ~ 0, fig. O we find the 



same sequences to be the branches of the separatrix which 
divides particular sequences of differentially rotating New- 
tonian bodies from one another. We mark this separatrix 
by a bold curve. The Maclaurin spheroids are character- 
ized by /3 = |, and the two branches discussed above can 
easily be identified. The Dyson ring sequence leads from a 
mass-shedding end f3 = across the Maclaurin sequence to 
the regime of toroidal bodies (at rp/r^ = 0). The remaining 
piece of the se paratri x is given by the sequence in fig. 4 of 
Ansorg et al. (|2003bh and leads also from the mass-shedding 
limit /3 — to the Maclaurin sequence. 

The separatrix composed of these segments looks quali- 
tatively similar to its polytropic counterpart in fig. [3] For 
differential rotation A > we find sequences of type C 
in the upper right and of type D in the lower left part of 
the panel, but neither type A nor type B. Note that the 
only type D curve displayed 'inside' the lower separatrix 
branches corresponds to yl = 0.1, just as the type C curve 
located opposite and next to the separation crossing point 
(rp/re = 0.17126,^= i). 

The curves printed in chain lines belong to the Maclau- 
rin sequence with rp/rc < 0.11160, P — as well as to the se- 
quence A^ in fig. 4 of Ansorg et al. l|2003bl '). i.e. to branches 
of separatrices dividing 'higher' solution classes from one 
another. 

In the remaining diagrams of fig. [5]we display sequences 
for A € {0,0.1,0.2, ...,1.5} and log Vax G {0.2,0.4,0.6}. 
As in fig. O the lowermost curves on the right part of the 
diagrams always correspond to rigid rotation A — 0, while 
the uppermost ones represent the configurations with the 
highest rotation parameter considered, A = 1.5. Again, bold 
lines describe the separatrix sequences, and the dotted lines 
in all diagrams of fig. [5] are sequences with (d^ Vd£'^)c = 0. 

For a given value log /imax the qualitative picture for ho- 
mogeneous bodies is similar to that for polytropes, i. e. we 
can identify sequences of types A and B for A < Aait 
and those of types C and D for A > Acrit. However, the 
onset of the spheroidal-toroidal transition as a function of 
log /imax is very different, see fig. [6] If one restricts oneself 
to differential rotation parameters A < AJ^jT^'^, then any se- 
quence of polytropic bodies starting at the spherical limit 
rp/rc ~ 1 ends in a mass-shedding limit /3 = with some 
value Tp/rc > (and consequently is of type A). This means 
that in this range of differential rotation we have a clear 
gap between the spheroidal and the toroidal class. The dis- 
tinction between the classes only disappears for A > A^"^^ 
in which case all four types of sequences can be found. In 
contrast, the solution space of differentially rotating homo- 
geneous bodies with A > does not allow at all a classi- 
fication with respect to spheroidal and toroidal classes. All 
four types of sequences can be identified in the small range 
< i < A^S^ = 0.29646, and for A > A^^ we find only 
type C and D sequences. That is to say that then any se- 
quence starting at the spherical limit rp/rc = 1 leads to the 
transition to the toroidal regime (and hence is of type C). 



5 SUMMARY AND DISCUSSION 

In this paper, we have constructed relativistic models of dif- 
ferentially rotating neutron stars with spheroidal topology 
(no hole) for broad ranges of degree of differential rota- 



Differentially Rotating Neutron Stars 7 




Figure 3. Sequences of differentially rotating polytropic fluid bodies (polytropic index N = 1) with prescribed maximal enthalpies 
log/iniax for A G {0, 0.1, 0.2, . . . , 1.5}. Uppermost and lowermost curves on the right part of the diagrams correspond to rotation 
parameters A = 1.5 and A = (rigid rotation) respectively. The solid bold curves represent the separatrix sequences with A = A^riti 
and the dotted curves defined by {d?h/dg'^)c = indicate that, for configurations to the left of these curves, the maximal enthalpy is 
assumed outside the coordinate origin. Note that the diagram for log hmax = displays Newtonian sequences. 
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Figure 4. For polytropic bodies (polytropic index A'^ = 1), the 
critical rotation parameter Acrit is plot against the maximal en- 
thalpy, log/imax- For parameters below this curve, each sequence 
starting at the spherical and static limit terminates at a mass- 
shedding limit (3 = 0. On the other hand, such sequences lead to 
the transition to the toroidal regime if the parameters chosen are 
located above the curve. 



tion A G {0, 0.1, 0.2, . . . , 1.5} and for maximal enthalpies 
log/iniax G {0,0.2,0.4,0.6}. We have considered homoge- 
neous and polytropic configurations (with polytropic index 
A'^ = 1) and adop ted the rotation law introduced by Ko- 
matsu et al. (|l989l '). 

From the results obtained we may conclude that the 
structure of disjoint classes in the solution space of rela- 
tivistic configurations is a special feature of the sub-space 
corresponding to rigid rotation. It turned out that this con- 
cept does not survive in the entire solution space of differ- 
entially rotating stars. Instead we find that from a critical 
rate of differential rotation on, a classification cannot be 
introduced. For differential rotation larger than this criti- 
cal rate, the previously disjoint classes have merged and we 
find continuous transitions from non-rotating and spherical 
to toroidal objects. 

The critical rotation rate depends strongly on the equa- 
tion of state being chosen. While the space of sufficiently 
weakly differentially rotating polytropes still permits the 
class structure in question, there is no non-vanishing rate 
of differential rotation for which such a classification makes 
sense in the case of homogeneous bodies. 

We have described properties of the general structure of 
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Figure 5. Sequences of difEerentially rotating homogeneous fluid bodies with prescribed maximal enthalpies log/imax for A £ 
{0, 0.1, 0.2, . . . , 1.5}. Uppermost and lowermost curves on the right part of the diagrams correspond to rotation parameters A = 1.5 
and A = (rigid rotation) respectively. As in fig. \3\ the solid bold curves represent the separatrix sequences with A = A^^n and the 
dotted curves are defined by (d^h/d(P)c = 0. Note that the diagram for loghmax = displays Newtonian sequences, some of which 
divide 'higher' solution classes from one another. See the text for more detail. 



the solution space in terms of four types of one-dimensional 
parameter sequences: 

• Sequences of type A exist for A < Aa-it and consist 
solely of spheroidal configurations. These sequences start at 
a static and spherical body and end at the mass shedding 
limit. Most of these configurations have the maximum en- 
thalpy located in the stellar center, but not all of them (for 
exceptions see e.g. fig. O panel for log /imax = 0.6, type A 
sequences in the vicinity of the separatrix curve). 

• Sequences of type B also exist for A < Aait- They start 
at the mass shedding limit and have a continuous parametric 
transition to a regime of bodies with toroidal topology. The 
branch with spheroidal topology configurations exists only 
for A larger than a threshold value. 

• Sequences of type C exist for A > ^crit- They start at 
a static and spherical body and possess a transition to ar- 
bitrarily thin rings. Consequently, type C sequences contain 
configurations with spheroidal and toroidal topology. 

• Sequences of type D also exist for A > Aait, but only 
in a narrow range of values of A. Their two extremities are 
configurations at the mass shedding limit. 

In the solution space, the type A- and B-sequences are 
clearly separated from the type C- and D-sequences by spe- 



cific separatrix curves which we determined in a rigorous 
numerical procedure. These separatrix curves, emphasized 
through bold solid lines in the figures (2] |3] and O have in 
most cases considered three branches that end in the regime 
of spheroidal bodies and one branch that leads to toroidal 
bodies. Note that in fig. [3] in the diagram for log /imax = 0.6, 
an exception can be seen, namely two branches leading to 
toroidal topology. 

The fact that the classification in classes fails in the 
realm of differentially rotating bodies means that a new in- 
terpretation of the criterion of the maximal mass for the 
distinction between a black hole and a rotati ng neutron star 
(see e.g. reasoning in Baumgart e et al . l|2000l )') becomes nec- 
essary. While Baumgarte et al. (|2000l ') discussed an increase 
of the maximum mass of differentially rotating stars, we find 
here that the family of sufficiently strongly differentially ro- 
tating bodies does not possess a maximum at all for the 
stars' (baryon and gravitational) masses. Evidently, one can 
always define one maximal value, restricting the calculation 
to simply connected bodies, as we did in this study, but that 
is not physically relevant since as the full family contains ar- 
bitrarily thin toroidal figures, the mass v alues may become 
arbitr arily large (see Ansorg et al. (|2003d ) and Fischer et al. 
(|2005l )'). In a forthcoming article, we shall come back to this 
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Figure 6. For homogeneous bodies, the critical rotation param- 
eter Acrit is plot against the maximal enthalpy, loghmax- As in 
fig. m parameters below or above this curve describe sequences 
starting at the spherical and static limit and leading to a mass- 
shedding limit or to the toroidal regime respectively. 



mass criterion and its application to differentially rotating 
fluids. 

In this article the complexity of the full solution space of 
differentially rotating fluid bodies was demonstrated, which 
could only be revealed by means of a precise and flexible 
code that can provide predictions which are relevant from 
the astrophysical point of view. Indeed, it should be noticed 
that the critical values obtained for A (--^ 0.7) are associated 
with a degree of differential rotation that is quite similar 
to what was obtained in dynamical simulations of merg- 
ing binar y neu tron stars based on such EOSs (see Shibata 
& Uryu ( 200C[ )'). This means that for the study of proper- 
ties of stationary configurations, considered as mimics of a 
specific phase before the collapse to a black hole, the exis- 
tence of multiple possible configurations has to be kept in 
mind during the numerical calculations, in order to deter- 
mine properly an actual sequence and not one built from 
configurations belonging to different types. 
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